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We apply the Chapman-Enskog procedure to derive hydrodynamic equations on 

<N. 

£jPy an arbitrary surface from the Boltzmann equation on the surface. 

< 

[>. ■ I. INTRODUCTION 

One of Boltzmann's greatest contributions to statistical mechanics is his kinetic equation 

o_| n 

^ ■ for gases The Chapman-Enskog (CE) procedure was first used to obtain the Navier- 

Stokes equations from the Boltzmann equation and later adapted to analyze the macro- 
dynamics of lattice gases 

a a. 

Recent work has extended the lattice gas automaton to 
simulate fluids on arbitrary surfaces jj], calling for methods of analyzing lattice gas methods 

>: 

■^J- ! defined on an arbitrary surface. A first step towards such methods is the extension of the CE 

■ expansion to a curved background space. In this paper we consider Boltzmann's equation 

in: 

on a curved surface, and use the CE expansion to obtain hydrodynamic equations on such 



surfaces. 

We consider the Boltzmann equation with the Bhatnagar, Gross and Krook (BGK) col- 



r 

lision operator |5|: 

5-H ' 



/ - / (0) 



d t f + v-Vf = - J - J —, (1) 

T 

where is the local equilibrium distribution, / is the single particle distribution function, 
and t is the typical time scale of the relaxation to local equilibrium. 

II. THE MOMENT EQUATIONS 

The fluid density, fluid velocity, and energy are defined in terms of / by averaging over 
all molecular velocities. 

p = J mfdv, pu = J vmfdv, p£ = J -m|v| 2 fdv. (2) 
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To hydro dynamic equations from ([Q) we multiply by the relevant moment and integrate over 
velocity. Using conservation of mass, momentum and energy in collisions gives: 

d tP + V • (pu) = 
<9t(pu) + V • P = 

d t (p£)+V -q = 0. (3) 

The pressure tensor P and heat flux vector q are defined by: 

P — J mvv/cfv, q = J -m\v\ 2 vfdv (4) 

The CE procedure enables one to derive constitutive equations giving P, and q in terms of 
p, u, and T for the fluid from successive approximations to Boltzmann's equation. 

III. THE CHAPMAN-ENSKOG EXPANSION 

The CE expansion is based on a physical picture of relaxation to equilibrium in which 
the fluid rapidly relaxes to a local equilibrium /(°) which depends upon space and time 
only through the variation of its moments. These moments then relax by hydrodynamic 
processes. The starting point is therefore an expansion of / around /^°^: 

/ = / (0) + + e¥ 2) + • • • (5) 

The terms are successive correction terms (for k > 0) and e is a formal parameter used 
to keep track of the order of the approximation. 

We assume that the macroscopic variables p, u, and T are given by the local equilibrium 
distribution f(°\ This implies: 

y"m/ (fc) rfv=0, f m/( fc Wv = 0, / f^\m\v\ 2 dv = 0, (6) 

for k > 1. Note that for the BGK equation this is equivalent to the requirement that the 
collision operator conserve mass, momentum and energy. 

The formal solution of the BGK equation may be written in terms of the linear differential 
operator T> = <9 t + v • V as / = (1 + rD) -1 /' '- We break up the operator T> into successive 
approximating operators D k : 

oo oo 
k=l k=l 
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Where the derivatives 9 t correspond to a hierarchy of time scales, with rapid relaxation to 



local equilibria on the shortest time scales, followed by hydrodynamic relaxation processes 
on longer timescales. 

Taylor expanding (1 + x)~ l we obtain a formal solution of the differential equation: 

k 

l + ^(_l)*y fc ( 



f 



k=i 



,i=i 



/ 



(0) 



(8) 



Equating this to the expansion of the distribution function ([5]) and collecting powers of 
e gives at first order p 1 ' = —rV^'f^. We compute the first order corrections to the 
equilibrium distribution and the viscous corrections to the momentum tensor and heat flux 
that they imply. 



IV. CONSTITUTIVE EQUATIONS 



The corrections to the pressure tensor and heat flux at kth order are defined as follows: 

P (k) = J mvv/ (fe W, q (fc) = J im|v| 2 v/ (fc) dv (9) 

To complete the derivation of the fluid equations we introduce an explicit form for the 
local equilibrium distribution function. We work in D dimensions to make the role of dimen- 
sionality explicit. We assume a Maxwell-Boltzmann equilibrium, although it is sufficient to 
assume that the local equilibrium distribution is Gallilean invariant and isotropic 
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The zeroth order energy density is: 
and the pressure tensor is: 



p(0) 



mvv 



pkT 



m 



1 + puu = pi + puu 



(10) 



(11) 



(12) 



where p= ^ is the scalar pressure and hence the gas has an ideal equation of state. The 
heat flux is: 



,(o) 



m v 



rf<®dv 



pu 



(D + 2)kT 1, ' 
2m + 2 |U| 



(13) 
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We will also need the heat flux tensor: 

Q(°) = [ mvvv/ (0) (fv = ■ u + puuu (14) 

J m 

where Qijki = Sij5ki + SikSji + 5u5jk is an isotropic fourth rank tensor. 

A. First Order: Euler's equations 

Defining the zeroth order terms for the pressure tensor P^ and the heat flux q(°) we 
obtain Euler's equations for an ideal fluid: 

d^p + V • (pu) = (15) 
d t (1) (pu) + V-P (0) =0 (16) 
d[ l) £ + V- (pq (0) ) = 0. (17) 

Simplifying and writing the momentum and energy equations in closed form in terms of the 
hydrodynamic fields gives: 

9 iU + u-Vu=-y, d t T + u-VT = TV -u (18) 

B. Second Order 

At second order we define and obtain the Navier-Stokes equations Q. 

pW = -rd tl J / (0) mvv(i o 
= _ r9fi p(0) _ rV . g(o) 

(19) 

eliminating the derivatives in favor of the hydrodynamic fields using the first order equations 
gives: 

pM = -I^{ V (u) + [V(u)f} (20) 

which is the correct form for the viscous stress in a Newtonian fluid with viscosity r] = 
rkTp/m. As the goal of the present paper is the derivation of incompressible fluid equations 
we do not compute the corresponding correction to the heat flux. 
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V. DIFFERENTIAL GEOMETRY ON TWO-DIMENSIONAL SURFACES 

In order to make this paper more self-contained, and to fix notation, we briefly re- 
view some relevant aspects of differential geometry. We consider an arbitrary surface, 
parametrized by coordinates a 1 for i — 1,2. Coordinates x l (a 1 ,a 2 ), i G {1,2,3} give the 
embedding of the surface in three-dimensional Euclidean space. The metric tensor on the 
surface, gij, is defined as: 

Edx 1 dx l , , 

MM' (21) 

i=i 

we define g to be the determinant of these components of the metric tensor, in terms of 
which we can express the Jacobian as J = ^fg. If e/ are the basis vectors for the surface 
coordinates u l the Christoffel symbols may be defined in terms of the derivatives of these 
basis vectors with respect to the surface coordinates: 

r M = e l ■ d e e k = r M = -g im {dig mk + d k g mi - d m g M ) (22) 

A contraction of the Christoffel symbol may be related to the determinant of the metric by 
T\ - = \d k log(g) and hence to the Jacobian JY\^ = d k J. 



VI. MOMENT EQUATIONS ON AN ARBITRARY SURFACE 

We consider the form of Boltzmann's equation on an arbitrary manifold. In the termi- 



nology of [S|, / is a relative scalar. To obtain an absolute scalar, in terms of which we may 
write an invariant form of the Boltzmann equation, we must multiply by the Jacobian of 
the coordinate system in which / is expressed. This reflects the physical fact that it is the 
number of particles in an elementary volume which is a scalar quantity, and if the volume 
form on a surface (or in an arbitrary coordinate system) varies from place to place, so must 
/, even for a constant density of particles. The invariant form of the Boltzmann equation is 
therefore: 

d t (Jf) + v i d l (fJ) = 0. (23) 

The macroscopic quantities p and e are also densities, and therefore relative scalars, while 
pu is a relative vector. Their definitions are unchanged from those given in equations (j2J). It 
should be noted that for any surface, the velocity vectors lie in a Euclidean tangent space and 
so we do not need to consider the variation of the volume form for the integrals over particle 
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velocity which define the macroscopic quantities. Finally, we will consider throughout a 
surface which does not change as a function of time and so d t J = and d t g l i = 0. 

A. The continuity equation 

The first moment equation is: 

J md t {Jf)dw + J v%{mfJ)dv = 0. (24) 

Expanding terms gives us: 

J d t p + Jditjw*) + diJ^pu 1 ) = 0. (25) 
The results of Section [V] allows us to cast equation (1251) as: 

Jd tP + Jdi(pu l ) + -lipinVy, = 0. (26) 

which is equivalent to: 

d tP +{pu% = U. (27) 

This is the continuity equation where the partial derivatives have been replaced by covariant 
derivatives on the surface. 

B. The momentum equation 

Multiplying by mv and integrating gives: 

J d t {J f)mv j dv + j v i v i d i (Jf)dv = 0. (28) 

We proceed as follows. First write: 

d t J Jfmv j dv = J d t (Jf)mv j dv+ j \jf)md t v j dv. (29) 

Usually the partial time derivative of the molecular velocity may be taken to be zero as the 
effects of particle interaction are taken into account only through the collision term. On 
an arbitrary manifold, in the absence of collisions or external forces the particles obey the 
geodesic equation of motion: 

d t v j + T j ik v i v k = (30) 



and we may write: 

dtj Jfmv j dv = J d t (Jf)mv j dv - j \j f)mT j ik v i v k dw (31) 
so that the first term in ( 128]) becomes: 

J d t (Jf)mv j dv = d t J Jfmv j dw + T j ik J {Jf)mv i v k dv. (32) 

For the spatial derivatives we proceed as for the continuity equation. We expand the partial 
derivative of J f: 

J v l v j di{fJ)dv = J v i v j Jd i fdv + J vVfdiJdv (33) 

once more, we may the results of Section |V] to write this in terms of the Christoffel symbols: 

J mv'viditffjdv = Jdi J mvWfdv + JT k kj J mvWfdv (34) 

Hence we may write the momentum equation on a fixed manifold which does not change as 
a function of time: 

d tP u j + diP ik + r j ik p ik + r k kj pv = o. (35) 

This is Cauchy's transport equation for momentum with the divergence of the momentum 
tensor replaced by the corresponding covariant derivative on the surface. 

C. The heat equation 

Multiplying Boltzmann's equation by |m|v| 2 = ]^mgjkV^v k and integrating gives: 

J ^mg jk v^v k d t (Jf)dv + J ^mg, jk v 1 v k v t d t {Jf)dv = 0. (36) 

again, assuming the surface is time independent so that d t J = and d t Qij = and using 
the same techniques as above for the last term here we obtain: 

Jd t (pS) + J 34 + JT^q* = (37) 

where we have defined the heat flux vector: 

q l = J l -mg jk v'v k v l fdw (38) 

Again we recognize this as the equation for the transport of heat with partial derivatives 
replaced by convariant derivatives. 

d t {p8) + q^ = (39) 
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VII. CHAPMAN-ENSKOG EXPANSION ON A MANIFOLD 

We proceed with the multiscale expansion of (E]), but now the expansion will be expressed 
in terms of derivatives on the manifold. The equilibrium distribution function is: 



./ 



(o) = y_ 



p f m \ D I 2 



m \2irkT 

the derivatives D^ k ' may now be written: 



Jexp 



-mgijiy 1 — u l ){yi — u^) 
2kf 



f)(k) = Q(k) + ^ + r jy 



(40) 



(41) 



As above, we only require the first correction to the equilibrium distribution, which may be 
written: 

/« = -t^W/P). (42) 
A. First order 

At first order we may compute the integrals which define the zeroth order pressure tensor, 
heat flux and heat flux tensor: 

p(o)« = 



q 



(0) 1 



mf^v' l v j dv = pg tj + pu l u j 

(p)ijk 



pu 



(D + 2)kT 1 ■ k 

h -Oi 7 M J M 

pkT 



Q(°) ijk = I mf ( - 0) v i v j v k dv = ^( u y* + u j g ik + u k g ij ) + pu i u j u k . 



m 



(43) 



Substituting these in the moment equations and simplifying for the case of an incompressible 
fluid with u l . A = gives the Eulers equations for the fluid: 



pdtu 1 + g ij djp + pddjU 1 + pu j u k Y) k = 
cfT 

— + u k d k T = Tu k k 



(44) 



and so the energy equation in the incompressible limit reduces to the statement that the 
convective derivative of temperature is zero. 

B. Second order 



At second order we must evaluate the correction to the stress tensor P^: 

P {1) = — t j vWfWmvVdv — — t J mv l v j (d? + v k d k + F^v 1 ) f {0) d^ 



(45) 
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once again we must correctly take into account geodesic motion when expressing P^ in 
terms of gradients of lower-order moments. In particular: 



mvWdPfWdv = d^P {0)13 - J /(°Wd t ( Vdv-y /<°W0 t ( Vdv 

= d\ X) pW j + r kl Q^ H + T j kl Q (0)ikl (46) 

so that: 

P^ 3 = -r{d^P^ 13 (47) 

by some lengthy but straightforward manipulations using our previously obtained first order 
quantities and their evolution equations we obtain: 

rpkT 



p(i) iJ 



(48) 



m 

which is recognizable as the generalization of the Newtonian viscous stress tensor where the 
usual partial derivatives have been replaced by covariant derivatives. 

VIII. CONCLUSION 

We have shown that the CE expansion on a manifold gives the Navier-Stokes equations in 
which the partial derivatives in the stress tensor are replaced by the covariant derivatives on 
the manifold. These are the correct fluid equations on the manifold, but these calculations 
offer us some additional physical insight. The fact that the covariant divergence of the 
pressure tensor has two terms involving the Christoffel symbols now has an explanation in 
terms of the underlying microscopic model. One of these terms arises because of variation 
of the volume form on the manifold. This term is also present in the divergence of the flux 
of scalar quantities. The second of these terms accounts for the fact that particles move on 
geodesies. This motivates a corresponding demonstration for the discrete fluid models on 
triangulations of 4| as follows: provided the triangulation approximates the volume form 
and geodesies of a manifold, the macrodynamics should approximate fluid dynamics on that 
manifold. We leave such a demonstration for future work. While the motivation for this 



work is two-dimensional manifolds relevant for [4] the extension to D dimensional manifolds 
is straightforward. This work is also of interest when considering fluids in the presence of 
background forces and interactions. Just as the Newtonian notion of force is replaced by 
the idea of motion on geodesies in General Relativity, one could include the effects of an 
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external force by a modification of the background geometry. Because the models defined 
in (4] also allow the geometry to change in response to the particle content the geodesies 
of the triangulation could perhaps encode the effects of the mean field forces present in an 
interacting fluid. 
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